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We present a simple phenomenological scaling theory for the pairing instability of a quantum 
critical metal. It can be viewed as a minimal generalization of the classical Bardeen-Cooper-Schrieffer 
theory of superconductivity for normal Fermi-liquid metals. We assume that attractive interactions 
are induced in the fermion system by an external 'bosonic glue' that is strongly retarded. Resting 
on the small Migdal parameter, all the required information from the fermion system needed to 
address the superconductivity enters through the pairing susceptibility. Asserting that the normal 
state is a strongly interacting quantum critical state of fermions, the form of this susceptibility is 
governed by conformal invariance and one only has the scaling dimension of the pair operator as 
free parameter. Within this scaling framework, conventional BCS theory appears as the 'marginal' 
case but it is now easily generalized to the (ir)relevant scaling regimes. In the relevant regime an 
algebraic singularity takes over from the BCS logarithm with the obvious effect that the pairing 
instability becomes stronger. However, it is more surprising that this effect is strongest for small 
couplings and small Migdal parameters, highlighting an unanticipated important role of retardation. 
Using exact forms for the finite temperature pair susceptibility from 1+lD conformal field theory as 
models, we study the transition temperatures, finding that the gap to transition temperature ratio's 
are generically large compared to the BCS case, showing however an opposite trend as function 
of the coupling strength compared to conventional Migdal-Eliashberg theory. We show that our 
scaling theory naturally produces the superconducting 'domes' surrounding the quantum critical 
points, even when the coupling to the glue itself is not changing at all. We argue that hidden 
relations will exist between the location of the cross-over lines to the Fermi-liquids away from the 
quantum critical points , and the detailed form of the dome when the glue strength is independent 
of the zero temperature control parameter. Finally, we discuss the behavior of the orbital limited 
upper critical magnetic field as function of the zero temperature coupling constant. Compared to 
the variation of the transition temperature, the critical field might show a much stronger variation 
pending the value of the dynamical critical exponent. 

PACS numbers: 



INTRODUCTION 



The 'mystery superconductors' of current interest share the property that their normal states are poorly understood 
'non Fermi- liquids'. Experiments reveal that these are governed by a scale invariance of their quantum dynamics. The 
best documented examples are found in the heavy fermion (HF) systems d i, i, H i, i, 0, i, Q . As function of pressure 
or magnetic field one can drive a magnetic phase transition to zero temperature. On both sides of this quantum critical 
point (QCP) one finds Fermi-liquids characterized by quasiparticle masses that tend to diverge at the QCP. At the 
QCP one finds a 'strange metal' revealing traits of scale invariance, while at a 'low' temperature a transition follows 
most often to a superconducting state with a maximum Tc right at the QCP. It is widely believed that a similar 
'fermionic quantum criticality' is governing the normal state in optimally doped cuprate high Tc superconductors. 
The best evidence is perhaps the 'Planckian' relaxation time observed in transport experiments ~ ^/(fcRT)[lol[lH 
indicating that this normal state has no knowledge of the scale Ep since in a Fermi-liquid r — {Ep/kBT)Tn- Very 
recently indications have been found that even the iron based superconductors might be governed by quantum critical 
normal states associated with a magnetic and/or structural zero temperature transition, giving rise to a novel scaling 
behavior of the electronic specific heat [l^ . 

The idea that superconductivity can be caused by a quantum phase transition involving a bosonic order parameter 
has a long history, starting with the marginal Fermi-liquid ideas of Varma 1J| in the context of cuprates of the late 



1980's and the ideas of spin-fluctuation driven heavy fermion superconductivity dating back to Lonzarich et 
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dealing with this subject that evolved since then 
departs from an assumption dating back to the seminal work of Herz in the 1970's[5l|. This involves the nature of the 
ultraviolet: at some relatively short time scale where the electron system has closely approached a Fermi-liquid the 
influence of the critical order parameter fluctuations become noticable. The Fermi surface and Fermi energy of this 
quasiparticle system can then be used as building blocks together with the bosonic field theory describing the critical 
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order parameter fluctuations to construct a perturbative framework dealing with the coupling between these fermionic- 
and bosonic sectors. The lowest order effect of this coupling is that the fermi gas of quasiparticles acts as a heat 
bath damping the bosonic order parameter fluctuations, with the effect that the effective space-time dimensionality of 
the bosonic field theory exceeds the upper critical dimension. These dressed order parameter fluctuations than 'back 
react' on the quasiparticle system causing 'singular' interactions in the Cooper channel, yielding in turn a rational for 
a generic 'high Tc' superconductivity at QCP's. 

The crucial assumption in this 'Herz philosophy' is that the fermion physics is eventually controlled by the Fermi 
gas. In the cases of empirical interest it is generally agreed that in the UV the interaction energies are much larger 
than the bare kinetic energies, while there is no obvious signature in the experiments for a renormalization flow 
that brings the system close to a weakly interacting fermion gas before entering the singular 'Herz' critical regime. 
From the theoretical side, the introduction of this UV Fermi gas can be viewed as an intuitive leap. The only truly 
fermionic state of matter that is understood mathematically is the Fermi gas and its perturbative 'derivative' (the 
Fermi liquid): the fermion sign problem makes it impossible to address fermionic matter in general mathematical 
terms }5 21]. However, very recently the 'grib of the Fermi-gas' has started to loosen specifically in the context of 
fermionic critical matter. A first step in this direction is the demonstration of proof of principle that truly critical 
fermionic states of matter can exist that have no knowledge whatever of the statistical Fermi energy scale: the 
fermionic Feynman backflow wavefunction Ansatzjssj. The substantive development is the recent work addressing 
fermion physics using the string theoretical AdS/CFT correspondence. It appears that this duality between quantum 
field theory and gravitational physics is capable of describing Fermi-liquids that emerge from a manifestly strongly 
interacting, critical ultraviolet (5J|. In another implementation, one finds an IR physics describing 'near' Fermi- liquids 
characterized by 'critical' Fermi surfaces [55l| controlled by an emergent conformal symmetry implying the absence of 
energy scales like the Fermi-energy 5^ 57|. 

This lengthy consideration is required to motivate the subject of this paper: a phenomenological scaling theory for 
a Bardeen-Cooper-Schrieffer (BCS) type superconductivity starting from the postulate that the normal state is not 
a Fermi-liquid, but instead a truly conformal fermionic state of matter. With 'BCS type' we mean the following: we 
assume as in BCS that besides the electron system a bosonic modes are present that cause attractive electron-electron 
interactions. This 'glue' is retarded in the sense that the characteristic energy scale of this external bosonic system ojb 
is small as compared to the ultraviolet cut-off scale of the quantum critical fermion system lUc ■ Having a small Migdal 
parameter, the glue-electron vertex corrections can then be ignored and the the effects of the glue are described in 
terms of the Migdal-jShashberg time dependent mean field theory, reducing to the static BCS mean field theory in the 
weak coupling limit [58]. All information coming from the electron system that is required for the pairing instability is 
encapsulated in the electronic pair susceptibility. Instead of using the Fermi gas pair susceptibility (as in conventional 
BCS), we rely on the fact that conformal invariance fixes the analytical form of this response function in terms of 
two free parameters: an overall UV cut-off scale (Tq) and the anomalous scaling dimension of the pair susceptibility, 
expressed in a dynamical critical exponent z and correlation function exponent rjp. The outcome is a scaling theory 
for superconductivity that is in essence very simple; much of the technical considerations that follow are dealing with 
details associated with modeling accurately the effects of the breaking of conformal invariance by temperature and 
the superconducting instability. This theory is however surprisingly economical in yielding phenomenological insights. 
Conventional BCS appears as a special 'marginal' case, and our main result is the generalized gap equation, Eq. (|10p . 
The surprise it reveals is the role of retardation: when the Migdal parameter lob/oJc is small (where the mathematical 
control is best) we find at small coupling constants A a completely different behavior compared to conventional BCS: 
the gap magnitude A becomes similar to the glue energy ujb- To illustrate the case with numbers, a moderate coupling 
to phonons like A = 0.3 with a frequency ojb = 50 meV will yield rather independently of scaling dimensions a gap of 
40 meV and a Tc of 100 Kelvin or so: these are numbers of relevance to cuprate superconductors! 

The theory has more in store. Incorporating the motive that on both sides of the quantum critical point heavy Fermi 
liquids emerge from the quantum critical metal as in the heavy fermion systems, we show that the superconducting 
'dome' surrounding the quantum critical point emerges naturally without changing the coupling to the bosonic glue. 
The form of this dome is governed by the correlation length, but we find via the pair susceptibility a direct relation 
with the effective mass of the quasiparticles of the Fermi-liquids. Last but not least, we analyze the orbital limiting 
upper critical magnetic field, finding out that pending the value of the dynamical critical exponent it can diverge very 
rapidly upon approaching the QCP, offering an explanation for the observations in the ferromagnetic URhGe heavy 
fermion superconductor 59]. 

The scaling phenomenology we present here is simple and obvious, but it appears to be overlooked so far. Earlier 
work by BalatskvjHoj. Sudbo61] and Yin and Chakravartv [6^ is similar in spirit but yet quite different. These 
authors depart from a Luttinger liquid type single particle propagators to compute the pair susceptibility from the 
bare fermion particle-particle loop. Although this leads to pair susceptibility similar (although not identical) to ours, 
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it is conceptually misleading since in any non Fermi-liquid, there is no such simple relation between two-point and 
four-point correlators. This is particularly well understood for conformal field theories: for the higher dimensional 
cases the AdS/CFT correspondence demonstrates that two point CFT correlators are determined by kinematics in 
AdS while the four- and higher point correlators require a tree level computation 63, 64, 65, 6^ 67, 68, 6^, 7^. More 



serious for the phenomenology, this older work ignores the role played by retardation; it is a-priori unclear whether 
one can construct a mathematically controlled scaling theory for BCS without the help of a small Migdal parameter. 

The remainder of this paper is organized as follows. In section II we review a somewhat unfamiliar formulation 
of the classic BCS theory that makes very explicit the role of the pair susceptibility. We then introduce the scaling 
forms for the pair susceptibilities as follow from conformal invariance. By crudely treating the modifications in the 
pair susceptibility at low energies associated with the presence of the pair condensate we obtain the new gap equation 
Eq. (fTU)) . This catches already the essence of the BCS superconductivity of quantum critical metals and we discuss its 
implications in detail. In section III we focus in on intricacies associated with determining the transition temperature. 
Conformal invariance is now broken and one needs to know the scaling functions in some detail. We use the exact 
results of 1-1-1 dimensional conformal field theory as a model to address these matters. In section IV we turn to the 
harder problem of modeling the crossover from the large energy critical pair susceptibility to the low energy, zero 
temperature infrared that is governed by conventional Bogoliubov fermions, as needed to devise a more accurate zero 
temperature gap equation. The casual reader might want to skip both sections. The moral is that information on 
the cross-over behavior of the pair susceptibility is required that is beyond simple scaling considerations to address 
what happens when the conformal invariance is broken either by temperature (as of relevance to the value of Tc) 
or by the presence of the BCS condensate (of relevance for the zero temperature gap). The conclusion will be that 
although the gross behaviors are not affected, it appears to be impossible to compute numbers like the gap to Tc 
ratio accurately since these are sensitive to the details of the cross-over behaviors. In section V we explore the theory 
away from the critical point, assuming that cross-overs follow to heavy Fermi-liquids, where we address the origin of 
the superconducting dome. Finally, in section VI we address the scaling behavior of the orbital limited upper critical 
field. 

BCS THEORY AND THE SCALING OF THE PAIR SUSCEPTIBILITY. 

Let us first revisit the backbone of Migdal-Eliashberg theory. We need a formulation that is avoiding the explicit 
references to the Fermi gas of the text book formulation, but it is of course well known how to accomplish this. Under 
the condition of strong retardation and small couplings, the effects of the glue are completely enumerated by the gap 
equationljH ignoring angular momentum channels {s,d waves, etcetera) for the time being, 

l-5X;ct(9 = 0,c^=.0,A,r) = 0, (1) 

where g is the effective coupling strength of the glue, while Xrct the zero frequency value of the real part of the 
retarded pair susceptibility at a temperature T in the presence of the gap A. This effective Xrot ^^^^ incorporates 
the effects of retardation. The textbooks with their focus on non-interacting electrons accomplish this in a rather 
indirect way, by putting constraints on momentum integrations. Retardation is however about time scales and the 
general way to incorporate retardation is by computing x'rct by employing the Kramers-Kronig relation starting from 
the imaginary part of the full electronic pair susceptibility Xp- For a glue characterized by a single frequency lob, 

XU^^0)^2V dc'M^. (2) 



with the fuh pair susceptibility given by the Kubo formula, 

Xp{q, u,) = -z ([fet(q, 0), 6(g, t)]) , (3) 



associated with the pair operator b'^{q,t) = ^l^)"^^ j^j^s |(^)- 

In the case of conventional superconductors the normal state is a Fermi-liquid, formed from (nearly) non-interacting 
quasiparticles. One can get away with a 'bare fermion loop' pair susceptibility. The specialty of this pair susceptibility 
is that its imaginary part is frequency independent at zero temperature. It extends up to the Fermi energy of the 
Fermi-liquid and from the unitary condition, 

/"OO 

x;{^)d^ = 1 (4) 
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FIG. 1: Illustration of the imaginary part of the pair susceptibility, comparing the relevant (Ising class), marginal (BCS case) 
and irrelevant (graphene class). The scaling exponent Ofp = respectively < cip < 1, cip = 0, Qp < 0. For the Ising 

class, the magnitude of the imaginary part of the pair susceptibility becomes larger and larger as one lowers the frequency. For 
the BCS case, the magnitude stays constant as the frequency is changed. For the graphene class, the magnitude decreases to 
zero in the low frequency infrared region. 



it follows that at zero temperature x"(a;) = A'o = 1/{2Ef)- In logarithmic accuracy the gap enters as the low 
frequency cut-off in Eq. ^ such that, 

x;.,(^ = 0,A,T = 0)= p-^ = -^log^, (5) 

and from Eq. ^ the famous BCS gap equation follows: A = 2uJBe^^^^, where A — g/Ep. 

This formulation of BCS has the benefit that it makes very explicit that all the information on the electron system 
required for the understanding of the pairing instability is encoded in the pair susceptibility. This is in turn a bosonic 
response function of the electron system since it involves the response of two fermions, much like the dynamical 
susceptibilities associated with charge- or spin densities. In addition one needs the fact that the pair density is a 
non-conserved quantity, in the same sense as a staggered magnetization. When the quantum system is conformal (i.e. 
the zero temperature quantum critical metal) the analytical form of the dynamical pair susceptibility is fixed at zero 
temperature by the requirement of invariance under scale transformations [t^ . 

xH^fulZ"[-[u^ + ^5fy-^^ (6) 

— 

as determined by the a-priori unknown unknown exponents "qp and z, the anomalous scaling dimension of the pair 
operator and the dynamical critical exponent, respectively. The normalization constant Z" is via the unitarity 
condition Eq.(|3]) determined by the UV cut-off scale lOc- Because we invoke a small Migdal parameter we are interested 
in the 'deep infrared' of the theory that is not very sensitive to the precise choice of this UV energy scale. A reasonable 
choice is the energy where the thermal de Broglie wavelength becomes of order of the electron separation, i.e. the 
Fermi energy of an equivalent system of non-interacting electrons. Defining Up = — and using Eq. (gl) with the 
cut-off scale we find, 

„ ^ I - ap 1 

sin(f ap) ^1-"^ ' ^ ^ 

observing that ctp < 1 in order for this function to be normalizable: this is the well known unitary bound on the 
operator dimensions. The real and imaginary parts of the zero temperature critical pair susceptibility are related by 
a phase angle f ap, 

^('^) ^ (cos(|ap) + i sin(|ap)) . (8) 

According to general conformal wisdoms, the pair operator is called irrelevant when ctp < such that %" increases 
with frequency, relevant when < < 1 when x" decreases with frequency and marginal when = 0, such that Xp 
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is frequency independent, see Fig 1. From this scaling perspective, the Fermi liquid pair operator is just the special 
marginal case, and the BCS superconductor with its logarithmically running coupling constant falls quite literally in 
the same category as the asymptotically free quantum chromo dynamics in 3+lD and the Kondo effect. Another 



familiar case is the pair susceptibility derived from the 'Dirac fermions' of graphene[73|, |7j| and transition metal 
dichalcogenides [tBI [To] characterized by ap = —1: in this 'irrelevant case' one needs a finite glue interaction to satisfy 
the instability criterium. 

The scaling behavior of the free fermion case is special and the pair operator in a general conformal fermionic state 
can be characterized by a scaling dimension that is any real number smaller than one. Obviously, the interesting case 
is the relevant one where ap > (Fig.l). Let us here consider the zero temperature gap equation. In Eq. ([6]) we have 
already fully specified Xp in the critical state. However, due to the zero temperature condensate the scale invariance 
is broken and the low frequency part of Xp will now be dominated by an emergent BCS spectrum including a s— or 
d— wave gap, Bogoliubov fermions and so forth. This will be discussed in detail in section V. Let us here introduce 
the gap in the BCS style by just assuming that the imaginary part of the pair susceptibility vanishes at energies less 
than A. Under this assumption the gap equation becomes, 

evaluating the integral this becomes our 'quantum critical gap equation' , 



with 



2A^^, (11) 



and A = g/ujc- The numerator (1 — ap) in A comes from the normalization constant Z", while the denominator ap 
from integrating over u. Notice that A has the same meaning of a conventional, say, dimensionless electron-phonon 
coupling constant. The dimensionful coupling constant g parametrizes the interaction strength between microscopic 
electrons and -lattice vibrations, and lUc has the same status as the Fermi-energy in a conventional metal as the energy 
scale that is required to balance g. We argued earlier that uic is of order of the bare Fermi energy and therefore it 
make sense to use here values for e.g. the electron-phonon coupling constant as quoted in the LDA literature. Notice 
however that for a given A the effective coupling constant A that appears in Eq. pop is decreasing when ap is becoming 
more relevant, i.e. when ap — > 1. From the frequency integral J dujuj~^-^^°'''\ one would anticipate that the gap would 
increase for a more relevant pair susceptibility. The unitary condition imposes however an extra condition on the 
pair susceptibility. These two compensating effects lead to the important result that the gap is rather sensitive to the 
relevancy of the pair susceptibility. All what really matters is whether the pair susceptibility is relevant rather than 
marginal or irrelevant, and the degree of the relevancy is remarkably unimportant. 

Eq. ljlOp is a quite different gap equation than the BCS one with its exponential dependence on the coupling A. The 
multiplicative structure associated with the Fermi-liquid is scaling wise quite special, while Eq. pU|) reflects directly 
the algebraic structure rooted in scale invariance. The surprise is that retardation acts quite differently when power 
laws are ruling. In Fig. (1) we show the dependence of the A/cub ratio on the coupling constant A, both for different 
Migdal parameters and fixed ofp, as well as for various scaling dimensions and the Migdal parameter fixed. The 
comparison with the BCS result shows that drastic changes happen already for small scaling dimensions ap especially 
in the small A regime. Our equation actually predicts that the gap to glue frequency ratio becomes of order one 
alrady for couplings that are as small as A = 0.1 when the Migdal parameter is small. To place this in the context of 
high Tc superconductivity, let us assume that the pairing glue in the cuprates is entirely rooted in the 'glue peak' at 



le peak at 



LOB ^ 50meV that is consistently detected photoemission, tunneling spectroscopy and optical spectroscopy 
. The electronic cut-off in the cuprates is likely of order Wc = 0.5 eV such that the Migdal parameter wb/wc — 0.1. A 
typical gap value is 40 meV and we read off Fig. 1 that we need A = 0.45 or 0.43 for ap = 3/4, 1/4 while using the BCS 
equation A = 1.1! Taking this serious implies that in principle one needs no more than a standard electron-phonon 
coupling to explain superconductivity at a high temperature in cuprate superconductors. Of course this does not 
solve the problem: although one gets a high Tc for free it still remains in the dark how to form a fermionic quantum 
critical state with a high cut-off energy, characterized by a relevant pair susceptibility. 

Eg.pi?]) is also very different from the gap equations obtained in the previous attempts to apply scaling theory 



to superconducting transition by Balatsky[60j, Sudbo|6l| and Yin and Chakravarty[62| . A crucial property of their 
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FIG. 2: (a)The ratio of gap to retardation frequency A/(2a;s) as a function of glue strength A for various retardation ranges 
ujb/i^c with fixed scaling dimension — 3/4. (b)The same plot, but with fixed retardation lob/i-Oc ~ 0.1 and various scaling 
dimensions ap. The dotted lines are the standard BCS result. 



results is that even in the relevant case one needs to exceed a critical value for A to find a superconducting instability. 
The present scaling theory is in this regard a more natural generalization of BCS theory, where the standard BCS is just 
the 'marginal end' of the relevant regime where the Cooper instability cannot be avoided for attractive interactions. 
The previous approaches 6^, 61 , start by considering the single particle spectral function, generalizing its analytic 
structure from simple poles to branch cuts. This way of thinking stems from the Fermi-liquid type assumption that 
the single particle Green's function is the only primary operator of the system, and all the higher point functions 
are secondary operators, to be determined by the single particle Green's function. But for critical systems, such 
assumptions are generally not to satisfied. It is well known for example from the AdS/CFT correspondence, that the 
four-point functions of strongly interacting conformal fields are much more complex than the combination of two-point 
functions ^1, |6j, |65|, |66j, IgtI-IgsI. [HqI . [70j . Our basic assumption is that the pair susceptibility is by itself a primary 
operator subjected to conformal invariance which is the most divergent operator at the critical point. 



DETERMINING THE TRANSITION TEMPERATURE. 



Let us now turn to finite temperatures. A complicating fact is that temperature breaks conformal invariance, 
since in the euclidean formulation of the field theory its effect is that the periodic imaginary time acquires a finite 
compactification radius = H/kBT. The pair susceptibility therefore acquires the finite size scaling form[73| 

Xiu;) EE x{q = 0,u;) = ZT'^^-^^^^^ (^) , (12) 

where $ is a universal scaling function and Z is a UV renormalization constant, while rjp and z are the anomalous 
scaling dimension of the pair operator and the dynamical critical exponent, respectively. At zero temperature this 
turns into the banch cut as shown in Eq.([6]), while in the opposite high temperature or hydrodynamical regime 
{huj « ksT) it takes the form[72| 

X{U0) = Z'T-(2-rfc)/^ 1 ^ (13) 

1 - ibJTrel 

where Trei ~ H/kBT. The crossover from the hydrodynamical- (Eq. to the high frequency coherent regime(Eq. 
I13p occurs at an energy ~ fc^T. The superconducting transition temperature is now determined by the gap equation 
through 1 — gx'y-cti^sTc) — 0. The problem is that x'ret is via the Kramers-Kronig transformation largely set by the 
cross-over regime in Xp- One needs the full solutions of the CFT's to determine the detailed form of $ in this crossover 
regime and these are not available in higher dimensions. 

In 1-l-lD these are however completely determined by conformal invariance, and for our present purposes these 
results might well represent a reasonable model since the gap equation is only sensitive to rather generic features of 
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the cross-over behavior. Given the exponents rjp and z, the exact result for the finite temperaturev" in 1+lD is well 
known and can be easily derived from the universal two-point correlator at an imaginary time r 72 1 

C{x,t) = Z——— — r~Tr^' (1^) 

(sin[7rj [t — ix)\ sin[7rj [t + ix)\) 

with 1 — 2s = — 27^- The analytical continuation to real time t it yields the real time two-point correlation function 

" ^ (^ sinh[^r(t - x)]i sinh[^r(t + x)])' ' ^^^^ 
with a Fourier transform corresponding to the dynamic structure factor 

nOO nOQ 

/ dx dtC{x,t)e''^''''-'^*\ (16) 



A convenient way to perform the Fourier transform is by factorizing C(x, t) into left-moving and right-moving modes, 
C{x, t) — C-{t — x)C+{t + x), to subsequently integrate over tztx. The result is 

^ ^ 1 „, Lo + k uj + k^„, u) — k bj ~ k , 
Sik, -) - Ze^J^^Bis + ^^,s - ^^)Bis + ^^,s - ^^), (17) 

where B is the beta function, and the overall numerical coefficient Z — 2^*^'^7r^(*^^^Z. The fiuctuation-dissipation 
theorem 

S{k,Lo) = j-^^x"{k,w) (18) 

then yields the imaginary part of the pair susceptibility, 

, sinh(7j^) , u; + k uj + k ^ , uj — k uj — k ^ , ^ 

X (k,uj) = Z ,,,^y, B(s + i——,s~i )B(s + i ,s-i , (19) 

The temperature and frequency dependencies of this function for fc = are illustrated in Fig. (3). Indeed x"(u^) ^ in 
a linear fashion with lo with a slope set by 1/T, while for u » T the temperature dependence drops out, recovering 
the power law. The crossover occurs at w ~ iksT /h where x" {lo) has a maximum. In the absence of retardation, the 
real part can be computed from the Kramers-Kronig transform, 

u)\ - — ( ~" s\n{2s-K-^)T(2s)r{2s-^) p (2„ . ^ j i^+Ii 2s- ^'l -I- H-i^ 1 - -^-11 

+ ik^ -k)) (20) 

where F is the generalized hypergeometric function. We did not manage to obtain an analytic form for the real part 
when retardation is included, and we use numerical results instead. 

When temperature goes to zero the limiting form of the beta function becomes, 

lim B{s + iu,s - iu) ^ ^^e-''''u^'-\ (21) 

U-+00 r(2s) 



and the imaginary part of the pair susceptibility Eq. (fT9l) acquires the power law form 



- ITL Z4-. (22) 



[r(2s)] 

Comparing this with Eq. ([7]) yields the normalization factor in terms of the cut-off scale 



(23) 

27r2(47r)"pwc 



Combining Eq.'s p]) . ((2|) . (fTQ]) . ((23|) . we obtain the equation determining the critical temperature 

\ LOc ) \2lob) V T, 



1_C'a(^V°7^V"V(^,.0. ,24, 
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FIG. 3: (a)Illustration of the imaginary part of the pair susceptibility ~x!' {k = 0, w,T) divided by the overall numerical factor 
Z, as a function of frequency uj for various temperatures. Hero we've chosen ap — 3/4, so s = 5/16. (b) The same plot, but x" 
is further divided by w. At zero temperature one has the power law scaling form. At finite temperature x"{^) goes to zero, 
as w goes to zero (x"(w)/w — > constant, as w — > 0), and approaches the same power law behavior at large frequency. As one 
increases temperature, the maximum of x"{'^) S^es down, and the corresponding uJmax shifts to larger frequency. 




FIG. 4: (a)The ratio of transition temperature to retardation frequency Tc/{2lub) as a function of glue strength A for various 
retardation ranges lob/cOc, with scaling dimension fixed Qp = 3/4. (b)The same plot, but fix the retardation cub /wc = 0.1 while 
varying the scaling dimensions Up. The dotted lines are the standard BCS result. The magnitude and dependence on glue 
strength and retardation are all similar to those of the gap. 



where 



Hy) = l --'^M2)(b(. + .^,.-^j^)) , (25) 

and X = uj/T. The overall coefficient is 

_ [r(2.)]^(l-a,) 



We plot in Fig. (4) the ratio of T,, to retardation frequency as function of glue strength, retardation and the scaling 
dimensions. One infers that the behavior of is very similar to that of the zero temperature gap, plotted in Fig. (2). 
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We observe that they are of the same order of magnitude Tc ~ A, and this can be understood from the behavior of 
x"/a; plotted in Fig. (3b). Since the large frequency behavior of x"(a;)/ti;'s are the same for different temperatures, all 
what matters is the low frequency part. The gap imposes a cut-off for the zero temperature x"(w)/a;, and its value 
is determined such that the area under this curve including the low frequency cut-off, is the same as the area under 
the curve corresponding to T^. without a cut-off: by inspecting Fig. (3b) one infers directly that the gap and Tc will be 
of the same order. The same logic is actually at work in the standard BCS case. The finite temperature Fermi gas 
susceptibility is x"{^) = tanh(i/3w) [ti'I, and the familiar Tc equation follows, 

1-A / — tanh(-/3u;) = 0, (27) 
Jo ^ 4 

such that Tc ~ l.liujBe^^^^ , of the same order as the BCS gap A = 2ti>Be^^/'^. Now the effect of temperature is 
encoded in the tanh function. Although the Fermi-gas is not truly conformal. It is easy to check that this 'fermionic' 
tanh factor adds a temperature dependence to the x" that is nearly indistinguishable from what one obtains from the 
truly conformal marginal case that one obtains by setting s = 1/2 in Eq. (fT9)) . 

We notice that conformal invariance imposes severe constraints on the finite temperature behavior of the pair 
susceptibility, thereby simplifying the calculation of T^. In the 1 -I- 1-dimensional 'model' nearly everything is fixed by 
conformal invariance. The only free parameters that enter the calculation are the scaling dimension ap, the cut-off 
scale ujc and the glue quantities. As we will now argue the situation is actually much less straightforward for the 
zero temperature gap because this involves a detailed knowledge of the crossover to the physics of the superconductor 
ruling the low energy realms. 



THE GAP EQUATION: GLUING A QUANTUM CRITICAL METAL TO A BCS SUPERCONDUCTOR. 

It is part of our postulate that when superconductivity sets in BCS 'normalcy' returns at low energies in the form 
of the sharp Bogoliubov fermions and so forth. Regardless the critical nature of the normal state, the scale invariance 
gets broken by the instability where the charge 2e Cooper pairs form, and this stable fixed point also dictates the 
nature of the low lying excitations. However, we are dealing with the same basic problem as in the previous section: 
in the absence of a solution to the full, unknown theory it is impossible to address the precise nature of the cross-over 
regime between the BCS scaling limit and the critical state at high energy. This information is however required to 
further improve the gap equation Eq. (jlOp of section II that was derived by crudely modeling x" in the presence of 
the superconducting condensate. 

So much is clear that the crossover scale itself is set by the gap magnitude A. However, assuming that this affair has 
dealings with e.g. optimally doped cuprate superconductors, we can rest on experimental information: in optimally 
doped cuprates at low temperatures the coherent Bogoliubov fermions persist as bound states all the way to the gap 
maximum. Up to these energies it is therefore reasonable to assume that Xp is determined by the bare fermion loops, 
and this regime has to be smoothly connected to the branch cut form of the Xp higher energies. This implies 
that the standard BCS gap singularities have to be incorporated in our zero temperature pair susceptibility. As a 
final requirement, the pair susceptibility has to stay normalized according to Eq. (|4]), which significantly limits the 
modelling freedom. 

Let us first consider the case of an isotropic s-wave gap singularity. The high frequency modes are still critical, and 
therefore the high frequency limit of the imaginary part of the pair susceptibility is determined by, 

hm x"(c^,A,T = 0) = ^. (28) 

In the presence of the superconducting condensate, the low energy modes below the gap have their energy raised 
above the gap, since we require x"{^ < A, A, T = 0) = 0. The spectral weight is conserved according to Eq. ([4]), and 
since we assumed that the Bogoliubov excitations of the BCS fixed point survive at energies of order of the gap we 
need to incorporate a BCS s-wave type power law divergence right above the gap in the imaginary part of the pair 
susceptibility. The simplest function satisfying these conditions is, 

X"(^,A)=^— ^=^_= e(.;-2A), (29) 

\ V^^^ - (2A2) j 
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FIG. 5: Illustration of the imaginary part of the pair susceptibility without a gap and in the presence of an s- and d- wave gap, 
for (a) the BCS case, and (b) the quantum critical case (here we've plotted using the parameter = 3/4). In the absence of 
gap, x" is a constant (for BCS) or has a simple power law behavior (for critical fermions). In the presence of a s-wave gap, the 
states below the gap are gapped out and there is a power law singularity right above the gap. When there is a d-wave gap, the 
low frequency part (way below the gap) is governed by a Dirac cone structure, thus a linear susceptibility, while near the gap 
a van Hove singularity is at work, leading to logarithmic divergences on both sides. The high frequency region for both s- and 
d-wave gap goes over to the case without a gap. 
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FIG. 6: The ratio of the gap to retardation A/(2a;s) as a function of the glue strength A, for various retardation ranges, with 
(a) a s-wave gap and (b) a d-wave gap. Here we've chosen Up = 3/4. The dotted lines are the standard BCS result. The 
dependence on glue strength and and retardation is similar but the magnitude of the gap is much enhanced compared to the 
previous treatment of gap as a simple IR cutoff. The d-wave case is enhanced even more. 

with ^ = (1 — ap)u)c^^~°'''^ (ycc Fig. 5b). We notice in passing that the BCS gap corresponds to the case Qp = 0, 



Xbcs('^,A) = 



CO 



-.e{(j-2A). 



'^Ef ^uP- - (2A)2 

The quantum critical gap equation for the s-wave superconductor now becomes, 



1 - 2(1 - ap)A ( 



2ujb 



A 

COB 



dx 



{x^ - l)(i+"p)/2 



= 0. 



(30) 



(31) 
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Turning to the d-wave case the gap equation becomes necessarily a bit more comphcated since we have to account 
for massless Bogolubiov fermions. At low frequencies lo << 2 A the pair susceptibility is now governed by free fermion 
loops and the Dirac-cone structure in the spectrum leads to a linear frequency dependence in the pair susceptibility, 
x"(a;) = Aiuj. Near the gap, a logarithmic divergence is expected due to the Van Hove singularity, and therefore 

X"(i^) = -Aalog for < 2A, while x"H = ^3 log for uj > 2A, with q, a cutoff. 

When the frequency is high compared to the gap scale, the pair susceptibility has the scaling form x"{^) — A^uj^'^'' . 
Matching these regimes at 2 A — wi and 2 A + ^2, with < wi < 2 A and < 1^2 < 2ti;s — 2 A, and assuming continuity 
of the pair susceptibility both below and above the gap (see Fig. 5b), we arrive at the gap equation for the d-wave 
case, 

1 wo A X ^ n''^^" dx , l + \/^Tl 



2g ^2AJo l-xg2/(2A) ^ -1 + V^TT 

2AJo l + a;9^/(2A) -1 + ^/a^ + 1 ap ^ 

This contains a number of free parameters that are partially constrained by the spectral weight conservation. This 
however does not suffice to determine the gap uniquely. In the following we will make further choice of the parameters, 
to plot the gap. We choose the scaling dimension ap — 3/4, and the cut-off in the logarithm to be of order the square 
root of the gap, say Qc/V^A — 3, the width of the logarithmic region to be 20 percent of the magnitude of the gap on 
both sides of the gap, that is wi/(2A) = cl;2/(2A) = 0.2, the coefficient of the high frequency part A4 — l/(4a;c^''), 
and further define ^i/g^ ^ ^^/^^ ^ „^ 5 ^ ja ^^j^g ^±^^,c ^ \og^±^^,d ee '"^^^o.32|•7.2I/^'''^^ ^^^^ the 
corresponding d-wave gap equation reads, 



1-0^ — — (0.8d + 7.2- / ---log- 



2 \ ujc J J c Jq 1 — 9a; —l + y/x + l 



1.2-4 f dx 1 + xAF + T 4, 3 ,A,3,, 

+ 9 / — —log— ^1== + - 1.2-3- =0, 33 

c Jq l + 9x -1 + V2; + 1 3 ujB 

We plot in Fig. (6) the behavior of the gap function in the s- and d-wave cases, to be compared with the outcomes 
Fig. (2) of the approach taken in section II where the gap simply entered as an IR cut-off scale, Eq. One can 
see that in both cases the magnitude of the gap is enhanced by treating the singularity more carefully, while in the 
d-wave case this enhancement is even more pronounced than in the s-wave case. These effects can be understood in 
terms of the redistribution of the spectral weight, since the low frequency part is enhanced by the factor 1/w in the 
Kramers-Kronig frequency integral. The dependence of the gap on the glue strength and retardation does however 
not change significantly compared to what we found in section II, which can be understood from the fact that the 
gap depends on the combination X{2lljb /'^c)'"'' ■ One also notices in Fig. (6) that the magnitude of the gap saturates 
already at small A for modest retardation. This is an artifact of the modeling. In real system the power law (s-wave) 



or logarithmic (d-wave) spectral singularities will be damped (see e.g. [80|, |8l|, |82|, |83| ) , and the endpoints at finite A 
in Fig. (6) will turn into smooth functions.. 

The gap to Tc ratio is expected to be a number order unity number. However, it is quite sensitive to the details 
of the crossover regime between the high frequency critical behavior and the low frequency superconducting behavior 
as of relevance to the zero temperature gap. Numerically evaluating Eq.'s (|24l31l33p we obtain gap to Tc ratio's 
as indicated in Fig. (7). Different from the Migdal-Eliasbergh case we find that these ratio's are rather strongly 
dependent on both the Migdal- and the coupling parameter, while the ratio becomes large for small coupling, in 
striking contrast with conventional strong coupling superconductivity. Invariably we find the ratio to be larger than 
the weak coupling BCS case, reflecting the strongly dissipative nature of quantum critical states at finite temperature 
that plays apparently a similar role as the 'pair-breaking' phonon heat bath in conventional superconductors. 



AWAY FROM THE CRITICAL POINTS: THE SUPERCONDUCTING DOME VERSUS T*. 



Our scaling theory yields a simple and natural explanation for the superconducting domes surrounding the QCP's. 
This is usually explained in the Moriya-Herz-Milles framework 51, 84, 85, 8^ that asserts that the critical fluctuations 
of the bosonic order parameter turn into glue with singular strength while the Fermi-liquid is still in some sense 
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FIG. 7: (a)The gap to Tc ratio 2A/T'c as a function of glue strength A for various retardation ranges uos/i^c with fixed scahng 
dimension ap — 3/4, for s-wave pairing. The dotted hne is the standard BCS result, where 2A/Tc — 3.5. (b) The same plot 
for d-wave pairing. The gap to Tc ratio decreases with increasing glue strength and retardation for both s- and d-wave gap. 
The ratios for different retardation ranges approach the same constant as A — > 0. 




FIG. 8: Illustration of the imaginary part of the pair susceptibility away from the critical point. For lo > T* , it has the critical 
scaling behavior, while for cu < T* , it retains the BCS form. T* is the cross-over scale. The effective mass m* is identified as 
the magnitude of the imaginary part of the pair susceptibility in the BCS region. The gap A acts as a low energy cut-off, and 
the retardation 2ujb as a high energy cut-off. When T* lies between A and 2uib, as is the case shown above, both the critical 
modes and Fermi liquid modes contribute. When A > T*, only the critical modes contribute. When 2a;s < T* , only the Fermi 
liquid modes contribute. 



surviving. We instead assert that the glue is some external agent (e.g., the phonons but not necessarily so) that is 
blind to the critical point, but the fermionic criticality boosts the SC instability at the QCP according to Eq. (|10p . 
By studying in detail the variation of the SC properties in the vicinity of the QCP it should be possible to test our 
hypothesis. The data set that is required is not available in the literature and let us present here a crude sketch 
of what can be done. In at least some heavy fermion systemsfs^] a rather sudden cross-over is found between the 
high temperature critical state and a low temperature heavy Fermi-liquid, at a temperature T* ~ \S — Sc\'^^ , with v 
behaving like a correlation length exponent ^ ^ |5 — ^d"" as function of the zero temperature tuning parameter 6. 
Moving away from the QPT this means for the SC instability that an increasingly larger part of the frequency interval 
of x" below u!B is governed by the Fermi- liquid 'flow' with the effect that Tc decreases. We can crudely model this by 
asserting that the imaginary part of the pair susceptibility acquires the critical form for lu > T* and the Fermi-liquid 
form for lo < T* , while we impose that it is continuous at a; = T* . This model has the implication that the magnitude 
of x" the Fermi-liquid regime is determined by T* and r/p and we find explicitly that A'o oc m* cx |(5 — 5c l"'^^^^''''^- 
We notice that this should not be taken literally, since this cross-over behavior can be a priori more complicated. In 
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FIG. 9: The ratio of Tc to retardation as a function of the distance away from criticality (a) for various scaling exponent a^'s 
with A = Q.QQ, uiB /tjJc = 0-1, vz = 3/2, (b) for various glue strength A's with cjs/wc = 0.1, i>z = 3/2, Op = 5/6. (c) for various 
retardation over cut-off o's/o'c's with A = 0.06, uz = 3/2, ap = 5/6. (d) for various inverse Griineisen exponent vz's with 
A = 0.06, lob/uJc = 0.1, ap = 5/6. 




FIG. 10: (a):Thc superconducting transition temperature Tc as a function of the distance from the critical point, for given 
crossover temperature T* and retardation lob. The parameters are chosen as z; = 2, = 1/3, r/p = 1, A = 0.05, ljb /(.^c ~ 0.1. 
(b):The same plot for a different set of parameters z — 3,rj = 0.5, f = 1/2, A = 0.05,a;B/a;c = 0.1. In-between the two points 
Sc ± 5, at which the transition temperature coincides with the cross-over temperature Tc{Sc ± 5) = r*((5c ± 5), the critical 
temperature remains constant. For T* > 2wb, Tc decays exponentially. The schematic behavior of the eflfective mass m* is 
also included. It diverges when approaching the critical point. 
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fact, from thermodynamic scaling it is known 8^, 8§| that m* ~ \5 — 5c\"'''^~'^^- Fig. (8) would imply that ap = 1 — d/z. 
This is not implied by scaling. 

Given these assumptions, the gap equation away from the quantum critical point becomes, 

1 - 2.9 £ ^X^cs(^) + ^Xc.it(^)) = (34) 

We are interested in the superconducting transition temperature, which has been shown in the previous section to be 
approximately the gap magnitude Tc ~ A. The imaginary part of the pair susceptibility in the critical region has still 
the power law form Xcnti^) ^ ^" siTi{apTr /2)uj~"'p , while in the BCS region it is a constant determined by continuity 
at w = T* and therefore Xbcs(^) = ^" sin(ap7r/2)(T*)-">' . 

Consequently we find in the regime Tc < T* < 2lub the solution for the gap equation. 



Tc = 2ujbx^^ exp 



— x"'-'^-''"^ - i( — )"''a;''(2-'''')^ 



(35) 



where x'^^ = T* /{2lub). For T* < Tc a plateau is found since only the critical modes contribute to the pairing, while 
for T* > 2ujb the BCS exponent takes over since only the (heavy) Fermi-liquid quasiparticles contribute having as a 
consequence, 

T,^2.,..J-i^)^:^)^ (36) 

The outcomes are illustrated in Fig. (9,10). One notices in all cases that the dome shapes are concave with a 
tendency for a flat 'maximum'. This is automatically implied by our starting assumptions. When Tc is larger than 
T* only the critical regime is 'felt' by the pairing instability and when this criterium is satisfied Tc does not vary, 
explaining the flat maximum. When Tc starts to drop below T* the superconductivity gets gradually depressed 
because the Fermi-liquid regime increasingly contributes. Eventually, far out in the 'wings', one would still have 
superconductivity but with transition temperatures that become exponentially small. The domes reflect just the 
enhancement of the pairing instability by the critical fermion liquid relative to the Fermi-liquid. 

The trends seen in Fig. 9 are easily understood. When the scaling dimension ap is increasing, i.e. the pair operator 
is becoming more relevant, the maximum Tc increases while not much happens with the width of the dome (Fig. 9a), 
for the simple reason that the critical metal becomes more and more unstable towards the superconductor. When 
the coupling strength A increases one flnds in addition that the dome gets broader (Fig. 9b) because the 'contrast' 
between Fermi-liquid and quantum critical BCS is becoming less, illustrating the surprise that especially weakly 
coupled quantum critical superconductors are much better than their traditional cousins. The same moral is found 
back when the Migdal parameter is varied (Fig. 9c), illustrating that at very strong retardation the differences are 
the greatest. Finally, in Fig. (9d) the evolution of the domes are illustrated when one changes the exponents relating 
T* to the reduced coupling constant. We find that the dome changes from a quite 'box like' appearance to a 'peak' 
pending the value of vz. The mechanism can be deduced from Fig. 10, comparing the situation that the quantum 
critical 'wedge' is concave (fig. 10a, vz < 1) with a convex wedge (fig. 10b, i^z > 1). Because T* is varying more 
slowly in he latter case with the reduced coupling constant, the quantum critical regime becomes effectively broader 
with the effect that the quantum critical BCS keeps control over a wider coupling constant range. The trends in Fig.'s 
(9, 10) are quite generic and it would be interesting to find out whether by systematical experimental effort these 
behaviors can be falsified or confirmed. 



SPATIAL DEPENDENCE OF THE PAIR SUSCEPTIBILITY: UPPER CRITICAL FIELD 

Another experimental observable that should be quite revealing with regard to scaling behavior is the orbital limiting 
upper critical field. The orbital limiting field is set by the condition that the magnetic length becomes of order of the 
coherence length, and the latter relates to the 'time like' Tc merely by the dynamical critical exponent z. In more 
detail, assuming a gap of the form [qOI , 



A(r) = Aoexp (^-^^ , 



(37) 
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the linearized gap equation in the presence of an orbital limiting magnetic field becomes 91| . 
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Xo(r,/3)exp(-^)r'^-idr, 



(38) 



where rt'^~^ is the volume of the d — 1-dimensional unit sphere, / is the magnetic length related to the field by 
H = 0o/(27r/^) where (f>o = hc/e, while Ko{r, /S) is the real space pair susceptibility, which is the Fourier transform of 
x'jil, 93J. For free fermions, the real space pair susceptibility is (see eg. [OJI), 



1 



1 



4/3sinh(|^)^ 



(39) 



with a power law behavior KQ{r, P) ~ r"'' at short distances or low temperatures where r < (3vp, and an expo- 
nential decay at large distances or high temperature. Let us consider critical fermions at T = 0, such that the pair 
susceptibility has the power law form x(ci^) ~ u^i^^'n)/^ ^ The momentum dependence can be determined by replac- 
ing ui by , such that x(^) ~ k^^'^^'^\ It follows that the real space pair susceptibility has the power law form 
KQ{r,T = 0) ^ J x(fc) exp(ifc • r)d'^k ^ j,-(ii-2+»)p)_ Associate with the retardation a short distance cutoff tq, and 
assume a scaling 2ujb/^c = {ro/ac)^^, where Oc is the lattice constant. The magnetic length acts as a long distance 
cutoff and therefore, 



1 



Ch 



dr, 



(40) 



with the normalization factor Ch — 2z{l — ap)n ^' 



SO that (l/fi'i-i) /^^ K,,it{ry'^dr ~ to give 



the right scale. The zero temperature upper critical field has then the same form as the one for Tc except for the 
occurrence of z, 



2^Hc2(0) 



1 + - 



2tjf 



and it follows, 



Poac 



2/z 



(41) 



(42) 



In the BCS case one has Fc2(O)/(S0ofcF) = {T^/Epf, with B ~ 3.26 for d = 3 [91. The moral is obvious: in 
Lorentz-invariant {z = 1) systems the relation between Hc2 and Tc is the same as for standard BCS, but when the 
normal state is governed by a universality class characterized by z > 1, Hc2{0) will be amplified for a given Tc relative 
to conventional superconductors because Tc/uJc,Tc/ Ep <C 1. 

Modeling the variation of Hc2 in the vicinity of the QPT as in the previous paragraph, where the critical modes 
govern the short distance and BCS type behavior is recovered at large distance, while converting the cross-over 
temperature to a length scale r* , by T* /ujc ~ {r* / Qc)^^ , we find that Hc2 is determined by the equation. 
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j,d—2+rip 



^dr 



^r'^'^dr, 

r* 



(43) 



with the matching condition Ch = (r*) ^+''pC,'j. We find that one just has to replace the first two dynamic exponent 
z's in Eq. (pS)) by 2 while an extra factor of 2 has to be added to the second term in the exponent. 



Hc2 — 



'I'nac ^ 2,y ( 



27r 



2/z 



exp 



u(2-rip) 



l/(2-7,p) 



A 



(44) 



In the region where only the Fermi-liquid quasiparticles contribute, the upper critical field has still an exponential 
form. 



He 



t>oac 
2tt 



2WF 



2/z 



exp 



-2( 



2u;b X' 



f(2-r)p) 



LUc 



(2-r])A 



(45) 
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FIG. 11: The upper critical field Hc2 over Hq = <j!>oaJ^/(27r) as a function of the distance away from criticality (a) for various 
scaling exponent ap's with A = 0.06, cjs/ljc = — 1/2, z — 3, (b) for various glue strength A's with ojs/tJc = 0.1, = 

1/2,2 = 3,ap = 5/6, (c) for various v's with A = 0.06, lub/uJc — 0.1, ap = 5/6, z = 3, (d) for various retardation ranges with 
A = 0.06, u = 1/2, 2 = 3, ftp = 5/6. 



The dependence of Hc2 on various parameters is shown in Fig. 11, and one infers that Hc2 behaves in ways very 
similar Tc (Fig. 10). The interesting part is illustrated in Fig. (12b) where we plot Hc2/Hq — Tc/uic as a function 
of the distance away from the critical point for different dynamical exponent z's, keeping all other quantities fixed, 
defining Hq = (j)Qa~'^ / {2tt) . One infers that when z > 2, Hc2/Ho — Tc/uJc increases rapidly when approaching the 
critical point. 

Using a 'ferro mag netic' dynamical exponent z = 3 and a Griineisen exponent l/vz — 2/3 inspired on recent 
experiments 9^ 9y| as well as theoretical considerations [^, |51, 55, 84, 8^, 8^ we obtain the results in Fig. (12a). 
Compared to Tc, Hc2 peaks much more strongly towards the QCP. This is in remarkable qualitative agreement with 
the recent results by Levy et al. on the behavior of the orbital limiting field in URhGe exhibiting a ferromagnetic 
QCPfHol, where the highest Tc is about 0.5 K 97 1, while the upper critical field exceeds 28 T. I t ha s also been 
observed in noncentrosymmetric heavy fermion superconductors CeRhSis [osl , |99| and CelrSia [l03,[lQl|, where the 
Pauli limiting effect is suppressed due to lack of inversion center of the crystal structures and the orbital limiting 
effect plays the main role of pair breaki ng. Near the quantum critical points, Hc2 can be as high as about 30 K, 



although the zero field Tc is of order IK [lol Il03l |. This class of experiments c an be understood in our framework as 
resulting from the change of the scaling relation between Hc2 and Tc. (See also [l04j for a tentative explanation from 
the customary Hertz-Millis-Moriya perspective.) 



CONCLUSIONS 



Perhaps the real significance of the above arguments is no more than to supply a cartoon, a metaphor to train 
the minds on thinking about pairing instabilities in non Fermi-liquids. This scaling theory has the merit of being 
mathematically controlled, given the starting assumptions of the 'retarded glue' and conformal invariance. The Migdal 
parameter plays an identical role as in conventional BCS theory to yield a full control over the gluc-fermion system 
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FIG. 12: (a) Illustration of the different behavior of Tc and upper critical field Hc2 as the quantum critical point is approached. 
Hc2 increases much faster than Tc. Thus for a small Tc one can still have a large upper critical field. Here we plotted using the 
parameters A = 0.05, ujb/l^c = 0.1, = 1/2, z = 3, 77 = —1. (b) The difference Hci/ Hq — Tc/uJc as a function of the distance away 
from the critical point for different dynamical exponent z's. Here Hq = (j)oa~'^ / (27r), A — 0.06, lobI^^c = 0.1, vz = 0.5, ap = 0.4. 
For z = 2, the difference is 0. For z = 3,4, the difference is positive and increases rapidly when approaching the critical point. 
For the case with z = 1, the difference is negative. 



dynamics, while we trade in Fermi-liquid principle for the even greater powers of scale invariance. The outcomes are 
gap and Tc equations where the standard BCS/Eliasberg equations show up as quite special cases associated with 
the marginality of the pair operators of the Fermi gas. The difficulty is of course to demonstrate that these starting 
assumptions have dealings with either nature itself and/or microscopic theories of electron systems where they should 
show up as emergent phenomena at low energy. However, the same objections apply to much of the current thinking 
regarding superconducting instabilities at quantum critical points with their implicit referral to a hidden Fermi gas. 
In such considerations there is an automatism to assume that eventually the superconductivity has to be governed by 
Eliashberg type equations. At the least, the present analysis indicates that such equations are not divine as long as 
the Fermi-liquid is not detected directly. Stronger, in line with the present analysis one might wish to conclude that 
superconducting instabilities will be generically more muscular in any non-Fermi-liquid. The Fermi-liquid is singular 
in the regard that its degrees of freedom are stored in the Fermi-sea, and this basic physics is responsible for the 
exponential smallness of the gap in terms of the coupling constant. This exponential smallness should be alien to any 
non Fermi-liquid. 

How about experiment? Scaling theories have a special status in physics because they guide the analysis of experi- 
mental data in terms of a minimal a-priori knowledge other than scale invariance. The present theory has potentially 
the capacity to produce high quality empirical tests in the form of scaling collapses. However, there is a great incon- 
venience: one has to be able to vary the glue coupling strength, retardation parameters and so forth, at will to test 
the scaling structure of the equations. These are parameters associated with the materials themselves, and one runs 
into the standard difficulty that it is impossible to vary these in a controlled manner. What remain are the rather 
indirect strategies discussed in the last two sections: find out whether hidden relations exist between the detailed 
shape of the superconducting and the crossover lines; are there scaling relations between Hc2 and Tc as discussed in 
the last section? We look forward to experimental groups taking up this challenge. 

There appears to be one way to interrogate our starting assumptions in a very direct way by experiment. Inspired 
by theoretical work by Ferrell |l05l | and Scalapino [l06| . Anderson and Goldman showed quite some t ime ago IiOtI 



that the dynamical pai r sus ceptibility can be measured directly using the AC Josephon efi^ect - see also 108, 109l |. fo 



a recent review see ref. |llO| . It would be interesting to find out whether this technique can be improved to measure 
the pair susceptibility over the large frequency range, 'high' temperatures and high resolution to find out whether 
it has the conformal shape. It appears to us that the quantum critical heavy fermion superconductors offer in this 
regard better opportunities than e.g. the cuprates given their intrinsically much smaller energy scales. 

In conclusion, exploiting the motives of retardation and conformal invariance we have devised a phenomenological 
scaling theory for superconductivity that generalizes the usual BCS theory to non Fermi-liquid quantum critical 
metals. The most important message of this simple construction is that it demonstrates the limitations of the usual 
Fermi-liquid BCS theory. The exponential smallness of the gap in the coupling is just refiecting the 'asymptotic 
freedom' of the Fermi-liquid, and this is of course a very special case within the landscape of scaling behaviors. 
Considering the case that the pair operator is relevant, we find instead an 'algebraic' gap equation revealing that at 



18 



weak couplings and strong retardation the rules change drastically: as long as the electronic UV cut-off and the glue 
energy are large, one can expect high T^s already for quite weak electron-phonon like couplings. If our hypothesis 
turns out to be correct, this solves the problem of superconductivity at a high temperature although it remains to be 
explained why quantum critical normal states can form with the required properties. It is however not straightforward 
to device a critical test for our hypothesis. The problem is the usual one that pair susceptibilities, A's or a^F's, and 
so forth cannot be measured directly and one has to rely on imprecise modelling. However, it appears to us that 
'quantum critical BCS superconductivity' works so differently from the Fermi-liquid case that it eventually should be 
possible to nail it down in the laboratory. We hope that the sketches in the above will form a source of inspiration 
for future work. 
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